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I. INTRODUCTION 

General relativity (GR) has passed all experimental 
tests in the Solar System and in binary pulsars with flying 
colors [l|. Current observations mostly probe the weak- 
field/slow-motion regime of the theory (with the excep¬ 
tion of binary pulsars, where the orbital motion is non- 
relativistic but the individual binary members are com¬ 
pact objects), and some of the most interesting strong- 
held predictions of GR are still elusive and difficult to 
verify. Observational and theoretical issues with Ein¬ 
stein’s theory - including the dark matter and dark en¬ 
ergy problems, the origin of curvature singularities and 
the quest for an ultraviolet completion of GR - have mo¬ 
tivated strong efforts to develop modihed theories of grav¬ 
ity which differ from GR in the infrared and ultraviolet 
regimes, while being consistent with the stringent obser¬ 
vational bounds at intermediate energies 0. The search 
for unambiguous signatures of modihcations of GR in 
the strong-gravity regime is a major goal of several re¬ 
search helds, including cosmology Q, “standard” elec¬ 
tromagnetic astronomy [^, and Earth- and space-based 
gravitational-wave astronomy [^, Q . 

In this work we consider a class of modihcations of 
GR known as Horndeski gravity Q. This is the most 
general scalar-tensor theory with a single scalar yielding 
second-order held equations for the metric and the scalar 
held (see e.g. [1, for tensor-multiscalar theories, and 
[i3,[iil for multiscalar versions of Horndeski gravity). All 
the terms present in the action of Horndeski gravity have 
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been shown to be originating from Galileons, i.e. scalar- 
tensor models having Galilean symmetry in flat space- 
time [l^. “Generalized Galileon” theories in curved 
space-time in any number of dimensions were studied 
in [H, and shown to be equivalent to Horndeski gravity 
in four dimensions in [l3| . Furthermore, Horndeski grav¬ 
ity can be shown to emerge from a Kaluza-Klein com- 
pactihcation of higher-dimensional Lovelock gravity (see 
e.gM for an introduction to this topic, and for a discus¬ 
sion of the relation between exact solutions in Lovelock 
and Horndeski gravity). 

The equations of motion of Horndeski gravity can be 
derived from the action 

S = J2 f , ( 1 ) 

where 

C2=G2, (2) 

£3 = -GaDcA , (3) 

£4 = G 4 i? + G4x [(□</-)'-<^^.] , (4) 

£5 = 

-^[(□^)3 + 20y . (5) 

Here is the metric tensor, g = det{gp,i,), and R and 
G^v are the Ricci scalar and the Einstein tensor asso¬ 
ciated with respectively. We have introduced the 
functions Gi = Gi{(j},X), which depend only on the 
scalar field (j) and its kinetic energy X = —d^(j)d^(l)/2, 
and we use units such that the reduced Planck mass 
rupj = (SttG)-^ = 1. For brevity we have also defined the 
shorthand notation (^^...^ = , 

and Gicj) = Horndeski theor¬ 

ies are an interesting phenomenological playground for 
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strong-field gravity because they include as special cases 
all dark energy and modified gravity models with a single 
scalar degree of freedom: 

(1) the GR limit corresponds to Gi{(j),X) = 1/2, with 

G2 = Gg = Gs = 0; 

(2) when the only nonzero term is = F((/) 

we recover a scalar-tensor theory with nonminimal 
coupling of the form F{cj))R, and therefore Brans- 
Dicke theory and f{R) gravity are special cases of 
Horndeski gravity; 


Einstein-dilaton-Gauss-Bonnet (EdGB) 
i.e. a theory with action 

gravity. 

S = l d'^x^(^^R + X + a^)Rl^'^ 

, (6) 

where -b Ra/g-y^R"^^^ is the 

Gauss-Bonnet invariant, corresponds to setting 

G 2 =V-b8^(^^V2(3-lnV) , 

(7) 

Gg ^ 

(8) 

G4=i+4^(2)^(2-lnV) , 

(9) 

Gs =-4^WlnV , 

(10) 

where Ra/s^yS R^^ are the Riemann and Ricci 


tensors, and we have defined = d^^jdcj)^ [T^ : 


(4) a theory with nonminimal derivative coupling of 
the form 

S = j d^x^[(:R + 2r]X + pG>^''(j)^(j),-2Ao] (11) 

(see e.g. [l 6 l - l^ for cosmological studies of this ^pe 
of action) corresponds to the following choiccQ of 
the coupling functions (2ll |: 

G 2 = -2Ao -b 2tjX , ( 12 ) 

G 4 = c + , (13) 

Gg = Gs = 0 , (14) 

where Aq, 77 , C ^ind (3 are constants; 

(5) the Lagrang ian £2 corresponds to the k-essence 
field (22 - 1^ (and therefore part of the literature 
uses a different notation, where G 2 is denoted by 

K); 

( 6 ) the covariant Galileon of Ref. [ 2 ^ is recovered by 
setting G 2 = —C 2 X, Gg = —csXjM^, G 4 = 
M|j /2 —C 4 X^/M® and G 5 = ScgAT^/M®, where the 
Ci (i = 2,..., 5) are constants and M is a constant 
with dimensions of mass. 


^ A coupling of the form <f >can also be obtained by setting 
G 5 = —<j) and integrating by parts. 


In this paper we are interested in black hole (BH) solu¬ 
tions in Horndeski gravity. As one of the most striking 
strong-held predictions of GR, BHs are ideal astrophys- 
ical laboratories to test gravity in the strong-held regime. 
Various authors explored nonrotating BH solutions in 
special classes of Horndeski gravity. Rinaldi [ 2 ^ stud¬ 
ied BH solutions in theories with a nonminimal derivat¬ 
ive coupling to the Einstein tensor of the form 
Minamitsuji (2^ and Anabalon et al. [ 2 ^ found more gen¬ 
eral solutions by adding a cosmological constant. Kobay- 
ashi and Tanahashi studied BH solutions in a sub¬ 
class of Horndeski theories that is both shift symmetric 
(i.e., symmetric under ^ -b c, with c a constant) and 
rehection symmetric (i.e., symmetric under <f) —>■ —</>). 
Under these assumptions, the only nonzero terms in the 
action are £2 and £4. Theories with nonminimal derivat¬ 
ive coupling are both shift and rehection symmetric, and 
therefore they are a subclass of the theories considered 
in Ref. [U. 

The nonrotating BH solutions found in the works lis¬ 
ted above either reduce to the Schwarzschild solution or 
are not asymptotically hat. This is a consequence of the 
no hair theorem by Hui and Nicolis [ 2 ^, which states 
that vacuum, static, spherically symmetric, asymptotic¬ 
ally hat BHs have no-hair in Horndeski theories with s hift 
symmetry. As pointed out by Sotiriou and Zhou [30|, l3l| 
the theorem actually has a loophole: asymptotically hat 
solutions can exist for theories of the EdGB type with 
= (j) (these theories are still shift symmetric, be¬ 
cause the Gauss-Bonnet combination is a topological in¬ 
variant). The solutions found in Refs. [^|m 1 are effect¬ 
ively special cases of the nonrotating EdGB BH solutions 
studied by Kanti et al. [s^ , that were subsequently gener¬ 
alized to slow rotation in Refs. [33l - l35l| and to rapid rota¬ 
tion in Refs. Other possibilities to violate the no¬ 

hair theorems include adding a time dependence to the 
scalar (but not to the metric), as in the solution proposed 
by Babichev and Charmousis [s^jOT considering biscalar 
extensions of Horndeski _gravity [ll|. Reference [s^ ex¬ 
tended the solutions in to the charged case, allowing 
for a coupling of the derivative of the scalar held to the 
energy-momentum tensor of the Maxwell held. Refer¬ 
ence (40| argued that a phase transition to charged hairy 
BHs can be realized through a nonminimal derivative 
coupling to the Einstein tensor; in this case, however, 
the equations of motion were solved perturbatively. 

The key question we address in this paper is the fol¬ 
lowing: does rotation produce interesting violations of 
the no-hair theorem at leading order in a slow-rotation 
expansion? In experimental terms, could we possibly 
observe violations of the no-hair theorem via frame¬ 
dragging experiments? The conclusion of our analysis 
is that frame-dragging corrections are exactly the same 
as in GR for all of the Horndeski BH solutions that we 
analyzed, with the (already known) exception of BH solu¬ 
tions in EdGB gravity [SOl - l^ . We do not expect this 
result to hold at second order in rotation, and this will 
be the topic of a follow-up study. 
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The plan of the paper is as follows. In Sec. m we 
present the equations of motion for slowly rotating BH 
space-times in Horndeski gravity, and we carry out some 
basic sanity checks (in particular, we check that GR and 
EdGB gravity are recovered in the appropriate limits). 
The field equations themselves are rather lengthy, and 
they are listed in Appendix O for the reader’s conveni¬ 
ence. In Sec. lIIIl we study slowly rotating BHs in theories 
with a nonminimal derivative coupling with the Einstein 
tensor, finding that frame-dragging corrections are ex¬ 
actly the same as in GR. In Sec. m we provide argu¬ 
ments (based on the work of Refs. (291 - 1^ ') to support 


this no-hair result. Finally, in Sec. |V] we present some 
conclusions and point out directions for future work. 


II. THE EQUATIONS OF MOTION 

The equations of motion that follow from the action 0 
can be written schematically as £ap = 0 (from variations 
of the metric) and = 0 (from variations of the scalar 
field), where 


£ ol ^ — 


9aP 


G 2 + G2xXal3 — 


GsxXa/sOej) -f -zdapG^^cf)^ — G3(o 


+ GoipG^ -l- GixXa/sR 


+ G/^fi^gaP — GiaP + 


GiXxXap — -GixQaP 


(□(()^ — (/)^^,) -I- 2n(/)G4X<('a/3 — 2 V(q,[G4X^/3)GI</>] 


-I- V^[GiX(t>^G\(j)]gap i^[GiX(t>(a4’^P)] — ^^l[GiX(|^^4’ap\ — ‘2G4X(I>P^4‘'^a + G{G^xXap 


~ 2^5gap) + ‘^G^cj)^ (^pGa)^ — X^[G54>(^aGp)^] ■ 


— X^lG^tjj^Gap] ■^{RG^^’aP — RapG^4>^^ 


+ □(G5<^a/3) + V„V^(G5<^^'^) - 2V^V(„[G5<();3)1 + 5a/3[V^V,(G5((.'^") - □(G5<^.")]} 

- -^{G^xxXap - -^gapG5x)[{^(ff + - -[G5x{0(j))^(j)ap - 2V(Q,[G5x(n</')" 

+ iffa/3V'^[G5x(n<(>)"</>^]} - - V"[G5X</>(a</>;.a<(>'^/3)] + ^ [G5X</>a</>^a<(>'^,3]} 

+ 2 {^5x(</>^y</>a/3 + ‘20(j)4>atJ.4>^p) - V(^ [G5X<(>a)[G5X(/>cr<('/ij/^^'"] 

— 2X^[G5x04>(l>{a4’P)ti\ + ^^[G5x04'(j)ij,(j)ap]} , 


(15) 


£p = G2 p + V„(G2x<(>“) - G3a“ - V„(G3x</-“n<^) ” □'^G3^ + G^pR + - <j>%)G^xp 

+ V“[G4XX<(>a(□</<' - (!>%)] + V‘^{G4xcI^c.R) + 2n(G4xn<(>) - 2V“V'5(G4 X</<c,3) + GspG^pr^ 

+ GfG^p - ^G5xp[{a(t>)^ + - S^l^Dcj)] + V“[G5x<^c<^^.G^"] - iv“{G5xx</>a[(n<^)" 

+ 2(/.3^ - 3<(.2,n<(.]} - in[G5x(n</>)"] - V“V^[G5X</>'^a<^^/3] + ^□(G5X<(>^.) + V“V'5(G5X<^a^n</>) . (16) 


Here we have defined Gia = XaGi, Xap = SXfSgap, 
and /i(a/ 2 / 3 ) = {haf 2 p + fipha)/^. These equations 
apparently contain higher derivatives, but they can be 
shown to be of second order using appropriate identities 
(cf. e.g. Appendix B of 0). 


To investigate the properties of slowly rotating BH 
solutions in Horndeski grav ity we follow the approach 
developed by Hartle [4l|, ^|, in which rotational correc¬ 
tions to the static, spherically symmetric background are 
introduced within a perturbative framework. At linear 
order in the BH angular velocity Q, the metric can be 


written in the form 

dr^ 

ds^ = —A{r)dt^ -\ --^-|-r^(d0^-|-sin^ 9dip^) — 2ijj(r)dtdip, 

B{r) 

(17) 

where the frame-dra-gging function w(r) is of order £1. 

Kobayashi et al. [^, 1^ carried out a fully relativistic 
analysis of linear perturbations around static, nonrotat¬ 
ing, spherically symmetric backgrounds. As a prelimin¬ 
ary step for this perturbative analysis, they derived the 
equations of motion for general static, spherically sym¬ 
metric vacuum space-times. Here we generalize these res¬ 
ults to the slowly rotating case, deriving the equations of 
motion for the metric component w(r). We also gener¬ 
alize the analysis of Refs, [i^, by allowing the scalar 
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field to depend on the radial and time coordinates, since 
a nontrivial time dependence of (j) allows for the existence 
of hairy BHs [H]. Following Refs. [111,133 , we assume the 
scalar field to have the functional form; 

<j) = = qt + ilj{r) . (18) 

Then the kinetic energy X is independent of t: 

where the prime means differentiation with respect 
to the radial coordinate r. Then in (fTKjl . Xap = 
—Tp''^5'^a5'^p/2 — 125*aS*p■ The tt and rr components 
of Eq. (fT^ yield two equations 


B. Einstein-dilaton-Gauss-Bonnet gravity 


EdGB gravity [HI corresponds to the choice of 
Eqs. (|71l- (fTU]) . If the coupling is linear in the field - i.e. 

= a(j) as in nanii, so that the theory is shift sym¬ 
metric - and g = 0 , we get 


{SaBcj}' — r)uj" -f 


240 

12acl>'B' -f 8aB(j)" -f- Bcj)' 

r 




— AaB ( j ) ^ + 7: —;-77 —4 


A 


A' B' 


A B 


uj' = 0 . (26) 


If instead we use an exponential coupling of the form 
^ and we set g = 0 , the frame-dragging equation 

becomes 


ftt = 0 , (20) 

err = Q, ( 21 ) 

and the scalar field equation of motion (ITHl) in the back¬ 
ground metric (El) is given by 

£0 = 0 , ( 22 ) 

where the explicit form of the left-hand sides of (I20II - (I22II 
is quite lengthy, and it can be found in Appendix For 

a static scalar field (g = 0), Eqs. (I20l) - (l22]) reproduce 

the results obtained in [4^ |4J] ; for reflection-symmetric 
theories, they reduce to the results of [ 2 ll| . 

For slowly rotating solutions at linear order in the BH 
angular velocity, the only nonvanishing component of the 
equations of motion yields a second-order ordinary differ¬ 
ential equation for the variable a;(r): 

= 0 . (23) 

Again, the explicit form of the left-hand side can be found 
in Appendix [Al 

Taken together, Eqs. dinii-iiisi) provide a full descrip¬ 
tion of vacuum space-times at linear order in rotation. 
We now consider two special cases as sanity checks of the 
equations of motion. 


A. General relativity 


As stated in the introduction, the Einstein-Hilbert 
Lagrangian of GR corresponds to setting G 4 = 1/2 and 
all the other functions equal to zero. In this case the 
equation of motion for the function u){r) simply reads 


UJ 


// 




(24) 


in agreement with the frame-dragging equation found by 
Hartle [4l| . If the nonrotating background is the Schwar- 
zschild solution this further simplifies to 



2re*4.) + 




— oj 


D/ 

2 </"r + 6(j)' + rcj)' [ 3—+ 2 </' 




in agreement with the result of Ref. 


III. NONMINIMAL DERIVATIVE COUPLING 
TO THE EINSTEIN TENSOR 

In this section we apply the formalism derived above to 
rotating solutions in a class of Horndeski theories charac¬ 
terized by a nonminimal derivative coupling with the Ein¬ 
stein tensor of the form ED- The theory defined by this 
action is invariant under both shift symmetry (</) —>■ (f>+c) 
and reflection symmetry (</> —>■ —(j)). Shift symmetry al¬ 
lows us to write the equation of motion for the scalar 
field (j) as a current conservation equation HH113 : 

= 0 . (28) 

In particular, for the action ED, the conservation equa¬ 
tion ( 1 ^ reduces to 

(W'^" - = 0 . (29) 

Moreover, following [^ we shall parametrize our solu¬ 
tions in terms of three auxiliary functions: 


1 

II 

<1 

(30) 

-2A/?77-b C»7 + Mo 
^ ’ 2 A /?2 

(31) 

g{X)=2{C,-pX) . 

(32) 


Using this parametrization, BH configurations within 
this theory can be easily obtained with the following pro¬ 
cedure. The tt component of the equations of motion, 
Eq. ([20]), leads to the equation 


UJ 


-w' = 0 . 


r 


(25) _^^|.[xc;(l-r2j-(A))] = 0, (33) 
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which can be integrated with the solution 

Xg‘^{X)[l-r^F[X)]=C , (34) 


where C is a constant. Equation (IMl) determines X{r) al¬ 
gebraically. Then the metric function A{r) can be found 
by solving Eq. (ED, which yields 


[vA)' 


1 - r^A 

^ 1 - r‘^T{X) ■ 


(35) 


Finally, the metric function B(r) can be found from 
Eq. (IMJ: 


B{r) = —[l-r^nX)\A{r). (36) 

With the choice o, the frame-dragging equation for 
uj{r) has a particularly simple form; 


Gu)" + uj' 


g.x'A- 



A' 

'a 



(37) 


As an extension of Ref. [2lj , we now consider nonrotat¬ 
ing BH solutions of Eqs. (El-iIMl) in different subcases 
and investigate the slow-rotation corrections predicted by 
Eq. dill) for each of these solutions. 

Case 1: T = 0. One possibility to satisfy Eq. (PI) 
is to impose X(Xf) = 0, where following the notation 
of [2lj we define Xf to be the value of X for which 
X{Xf) = 0, and C = XfG'^^Xf). In this case, the metric 
components and the scalar field read 


A{r) 

B{r) 


r 


2Xt 


2^/ ^ 
=^A{r) , 
q 

q^-2XjA{r) 

A{r)B{r) 


V 2 
—r 

3/3 


(38) 

(39) 

(40) 


where /i is an integration constant. 

With a rescaling of the time variable q^ = 2Xf, 
Eqs. (Pl)-(|ini) represent a BH solution with an effective 
cosmological constant A = —rj/(5 and a nontrivial profile 
for the scalar field. Replacing this solution into Eq. (EZl) 
we find that uj{r) satisfies the same equation (I25() as in 
GR. The standard solution of this equation is 


W = Cl -I- ^ 


(41) 


where ci and C 2 are integration constants which can be 
fixed by imposing appropriate boundary conditions. 

Case 2: 0 = 0. Another class of solutions of Eq. (1341) 
corresponds to choosing G{Xg) = 0. In this case, from 
Eq. (l37l) we see that the coefficients of both ui" and w' 
vanish, and there are no corrections at linear order. 

Case 3: q = 0. Finally, we consider the case in which 
the scalar field is time independent (<7 = 0). Integration 
of the equations of motion for A{r) and B{r) leads to 


A{r) 


B{r) 


+ /3Ao (3/3 - Tyr^)] - 

V?(/3Ao + arctan 
+ ’ 

(2/3C - /3Aor2 -|- 


(42) 

(43) 


where again ^ is an integration constant, while for the 
scalar field we obtain: 




(/3Ao -b C>?) [r^iCv - /3Ao) -b 2/3Cr] ^ 
(/3 -I- Tyr^)^ A{r) 


Replacing the former expressions into Eq. dSil), we find 
that the frame-dragging function ui{r) satisfies once again 
the same equation (1^ as in GR. 


IV. WHY THE BALDNESS? 


The no-hair theorems for static, spherically symmetric 
BHs proved in Refs. [2^-[3l| rely crucially on shift sym¬ 
metry, which allows us to write the equation of motion 
for (j) as the conservation equation (1281) . In this section 
we discuss how these theorems can be generalized to the 
case where we consider first-order rotational corrections 
and time-dependent scalar fields of the form (1181) . In this 
case, we can show that the nontrivial components of 
are given by 
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r= Bii)' 


— G 2 X 


BiP' 


A' 

T 


G 3 X 


2B 

r 




2r2 A 


A 


BA' 

-^^G3X 


A 2r2 


A' 

T 

■,2 


B- 1 


2j,/2 


-Bijj' 


Br 
'2B A 


G 4 X — 


2B^^ 

r 

?2, 


A' 

'A 


G 4 XX 


B^^'A' 


(45) 


q 


G 2 X — 


„ B fB' A 4 , , 


G3X-- B' 
r 


B- 1 


G 4 X 


2B^'il)' 


B' A 1 , ,, 
s'+ T+r 


G 4 XX + ^ 




[2r 

G5X 


B^'tP 


12 


2r2 

r 

’ 


2r+(^ + f )^' 


2A 

r 


Bib’ 

G4X - ^Gsx 
2r 


Gsxx 


+ I 2 G 4 X + 


B^' 


Gk- 


2r'!/)' 


S 

]4 


2AiP' 


rG'-, 


g; 


• (46) 


For shift-symmetric theories, Gi = Gi{X). These ex¬ 
pressions can be used to extend the no-hair theorems of 
Refs. [2^-[3l| to the cases considered in this paper. 

For clarity and completeness, let us begin with a short 
summary of the original proof given in (with the 
amendments of Refs. l^b 


which can be easily integrated. The solution is 
= K, where K is an integration constant. At 
the horizon the areal radius r cannot be zero. This 
implies that K = 0, and therefore that 

r = 0 V r . (49) 


A. A review of the no-hair theorem for nonrotating 

black holes with a time-independent scalar field , rr^i j. rr u u j.- 11 j-j. 

4. ihe current J can be schematically written as 


The no-hair theorem of Ref. [ 2 ^ applies to static, 
spherically symmetric, asymptotically flat solutions in 
shift-symmetric theories. It consists of the following line 
of reasoning: 

1. Assuming that the scalar field '!/'(''’) has the same 
symmetries as the metric (the time-dependent 
scalar field of [s^ obviously violates this first as¬ 
sumption) , the only nonvanishing component of 
for a spherically symmetric background is J*", i.e. 
JM = (J*-, 0,0,0). 


2. Given a spherically symmetric space-time, defined 
by the line element (1171) with uj{r) = 0, we require 
to remain finite at the horizon rh- Since 


J2 


B 


(47) 


and i? —> 0 for r ^ rh, this regularity condition 
implies that J*" = 0 at the horizon. 


3. For a spherically symmetric space-time, the conser¬ 
vation equation (1^ reduces to 

-^d^{./^gj'') = drr+ -r =0 , (48) 


r = B^'F{g,g',g",i;') , (50) 


where F" is a generic function of the metric, its hrst 
and second derivatives, and tp'. At spatial infinity, 
asymptotic flatness implies that R —>■ 1 and ip' —>■ 0, 
while F tends to a nonzero constant. This last 
condition is dictated by the requirement that the 
scalar field’s kinetic energy should have the stand¬ 
ard form: in the weak-field limit, the action con¬ 
tains a term that is quadratic in the field derivatives 
and —>■ dfi<p, up to an overall constant of normal¬ 
ization. If we now move “inward” towards the ho¬ 
rizon, by continuity F and B will still be nonzero, 
and therefore J” ^ 0, which contradicts Eq. (H^ . 
This contradiction can be avoided if '0' = 0 for any 
choice of r, which fixes ip = constant or (without 
loss of generality, since the theory is shift symmet¬ 
ric) 0 = 0. 


Sotiriou and Zhou [s^, Ej pointed out a loophole in 
the last step of this proof. For Horndeski gravity the¬ 
ories with shift symmetry, the conserved current can be 
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written as 


r 


BG2xi)' + 


2 


f A' 1 
r \ A Br 



Gax 


/A' 1\ 

r \ ^ r J 
A' /3B- 1 
2r2 ~A V~B 


G4XX 
^ Gsx + 


A' 54^/4 
T 2r2 


Gsxx • 


(51) 


Depending on the particular form of the coupling func¬ 
tions Gi we have essentially two options 0 

(a) depends linearly on ■0'. This is the case con¬ 
sidered in Ref. [^, for which F —G 2 X as 
r —>■ 00 . 


(b) contains terms which are independent of ip', but 
no negative powers of ip'. 

This second case represents a loophole for the no-hair 
theorem of Ref. [ 2 ^. Indeed, in this case the asymptotic 
behavior of F is not trivially determined. 

This is illustrated most clearly by looking at two spe¬ 
cific examples: EdGB gravity and theories with nonmin- 
imal derivative coupling to the Einstein tensor. 

In the first case the conserved current reduces to 

4dGB = -BiP’ - , (52) 

where we specialized to a linear coupling function ^ [p) = 
ap in Eq. ©, so that the theory becomes shift sym¬ 
metric (recall that the Gauss-Bonnet combination is a 
topological invariant). The current (j52p contains a term 
independent of p' as in case (b) above, corresponding to 
the loophole pointed out in Refs. [10, [3l|. The current 
vanishes at infinity, but for smaller radii the choice of F 
is nontrivial and leads to scalar hair growth. 

For the nonminimal derivative coupling theory we have 
instead 


= BP' 


^ 2R 

2r]-\ - 

r 


y A Br r ) 


P 


(53) 


B. Extension to slow-rotation and time-dependent 
scalar fields 

What is crucial for the present work is that the argu¬ 
ments above apply also to rotating BH solutions at linear 
order in rotation. This is because, as argued in Ref. [30l| . 
the scalar field p (like all scalar quantities) is affected 
by rotation only at second order, and therefore the ex¬ 
pression ()45l) for the current ,F remains unchanged at 
linear order. Similarly, J® is still equal to zero at linear 
order. The component J'^ acquires a nonzero value pro¬ 
portional to the BH angular momentum; however J‘^ is 
independent of tp, and therefore it does not contribute to 
the current conservation equation ( 1 ^ . 

At first sight, the fact that no-hair theorems still hold 
true at linear order in rotation even for time-dependent 
scalar fields may be surprising. However this no-hair 
property can be proved through a simple extension of 
the arguments valid for static, nonrotating solutions. Let 
us extend the original argument to theories with time- 
dependent scalar fields of the form (fTSl) : 

1. When p has the form (11811 the current has a nonzero 
time component, i.e. = (J’', 0, 0, J*), and its 
norm becomes 

j2 ^ ^ (54) 

B 

2 . By imposing regularity at the horizon, where A —>• 
0 , H —>■ 0 , we conclude that J’’ —0 as r — Vh- 
This is true as long as J* does not diverge in the 
limit r —>■ r^, i.e., as long as the quantity in square 
brackets in the last line of Eq. (HSl) is finite. For 
reflection-symmetric theories (G 3 = G 5 = 0), this 
latter requirement simplifies to the condition that 
{B/A)' should be finite [llj. 

3. In principle, the current conservation equation (j48p 
acquires an extra term because J* 7 ^ 0: 

drr+ -r+ dtf = G . (55) 

r 


This expression for the current falls into case (a) above. 
The current depends linearly onp', F ^ —2r] for r —>• 00 , 
and F stays finite even at finite radii by continuity, as 
required by the arguments of [l^, so we are forced to 
set p' = 0 and p is a. constant, which can be set to 
zero. Asymptotic flatness was of course a key ingredi¬ 
ent in these arguments. Hairy solutions in theories with 
nonminimal derivative coupling are not asymptotically 
flat (see e.g. [ 261 - 1 ^ '). 


^ A third case where J’’ contains negative powers of iji' can be 
excluded because it generally corresponds to theories that would 
not admit flat space with a trivial scalar conflguration as a solu¬ 
tion, leading to violations of local Lorentz symmetry l3ll . 


However Eq. (14611 shows that in the present case 
J* is independent of time, so this term vanishes: 
dtJ*" = 0. Following the reasoning below Eq. (H51l . 
we conclude that J’" = 0 for all r even for scalar 
fields with a linear time dependence. Note that 
for a time-dependent scalar field, in general, the tr 
component of the gravitational equations £tr = 0 
may be nontrivial, indicating the existence of an en¬ 
ergy flux in the radial direction. However Ref. 
showed that, for the linear-in-time ansatz m, 
Str is proportional to J*" under the assumptions 
of diffeomorphism invariance and shift symmetry. 
Therefore the condition ,F = {) always ensures that 
£tr = 0 : the linear time dependence (11811 does not 
give rise to an energy flux in the radial direction. 













4. The current (1451) has the form dsni), where 
F{g,g',g "is an unspecified function. This al¬ 
lows us to borrow in its entirety the reasoning of 
Ref. [m. We can exclude cases where J’' contains 
negative powers of ip'. When all terms in J'' contain 
positive powers oi ip', ip' = 0 for all r and the no¬ 
hair theorem of applies. The only exception is 
the case where J*" contains one or more terms with 
no dependence on ip', but no terms with negative 
powers of ip'] and then, following Sec. IIB of 
shift symmetry and Lovelock’s theorem imply that 
the action must contain a term proportional to the 
Gauss-Bonnet invariant. 


This generalized no-hair theorem can be used to justify 
the absence of corrections to GR at linear order that we 
found in Sec. m For a theory with nonminimal derivative 
coupling to the Einstein tensor, the nonzero components 
of the current can be obtained by specializing Eqs. dUD- 
(l46l) . with the result 


rG, = Bip' 





A 


Jgs — 277-1- 


2B fl- B 


Br 




(56) 

(57) 


The J*" component is identical to the static case of 
Eq. (15^ . it does not contain any 7/;'-independent terms, 
and the no-hair theorem of [ 2 ^ implies that asymptotic¬ 
ally flat solutions must be the same as GR. 

In conclusion, the only no-hair violations at linear or¬ 
der in rotation when the scalar field depends linearly on 
time and when we require asymptotic flatness can occur 
in one of two cases: 


(i) if the scalar field has a linear coupling to the Gauss- 
Bonnet invariant, or 

(ii) if, as proposed in Ref. [s^ , the field equations of the 
theory guarantee that the current vanishes identic¬ 
ally (J’' = 0) because F{g,g',g",ip') = 0 as a con¬ 
sequence of the field equations. Note that this is 
only possible for special forms of the functions Gi, 
and that the scalar field must then be time depend¬ 
ent (i.e., it must violate some of the symmetries of 
the metric) in order to be regular at the horizon. 


V. CONCLUSIONS 

In this work we studied leading-order rotational cor¬ 
rections to a broad class of BH solutions in Horndeski 

S ivity. With the known exception of EdGB gravity [ 13 - 
, we have found that the frame-dragging function U!{r), 
which describes the leading-order rotational corrections, 
is exactly the same as in GR for all of the Horndeski BH 
solutions known in the literature. This result applies even 
to asymptotically flat solutions that violate the no-hair 


theorems by requiring the scalar field to be time depend¬ 
ent (so that the scalar field does not respect the same 
symmetries as the metric), as proposed in Ref. 

The formalism developed in this paper can be extended 
in various directions. Eirst of all, the no-hair theorem 
proved in Sec. m at first order in rotation is not expected 
to hold at second order, where the continuity equation 
will be modified. Calculations of BH solutions at second 
order in rotation, along the lines of [U, [s^l , are already 
underway (45l| . 

Even for nonrotating Horndeski BHs, studies of sta¬ 
bility and perturbative dynamics (as encoded in their 
quasinormal mode spectrum: see e.g. 14611 for a review) 
are still in their infancy. One of us |47l | studied mass¬ 
less scalar field perturbations of static BH solutions in 
theories with field derivative coupling to the Einstein 
tensor. More in general, gravitational perturbations of 
static, nonrotating space-times can be explored using the 
formalism developed in Refs. [43,I3- The present work 
lays the foundations to study quasinormal modes and 
look for super-radiant instabilities using the slow-rotation 
perturbative techniques reviewed, e.g., in Ref. [4^. 

Another important extension concerns compact stars 
in Horndeski gravity. Slowly rotating compact stars in 
EdGB gravity were studied in [i^l ■ Cisterna et al. [131 in¬ 
vestigated compact objects in theories with a nonminimal 
derivative coupling of the scalar field with the Einstein 
tensor. Our formalism can be extended relatively easily 
to study compact stars in broader classes of Horndeski 
gravity, and to understand whether genuine strong-held 
deviations from GR (similar to the “spontaneous scalar- 
ization” phenomena proposed by Damour and Esposito- 
Earese [51|) can occur in some sectors of the Horndeski 
gravity action, see e.g. [13] for recent work in this direc¬ 
tion. 
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Appendix A: Field equations 


In this appendix we list the left-hand side of the field 
equations. For clarity, we split all of the left-hand sides 
of the field equations as a sum of two contributions, so 
that the case of time-independent scalar fields can more 

I 


easily be recovered by setting q = 0: 

+ (Al) 

= + (A2) 

Let us remark that the equations of motion still depend 
on the specific form of the G^’s, which are functions of 
the kinetic energy m, and therefore may contain q- 
dependent terms; therefore we must evaluate all of the 
functions Gi at g = 0 to recover the time-independent 
limit. The explicit forms of the various terms are 




W + 2 G,x - ^ f ^ + fid a - ^ f ^ + 4^' + 2^ I 


B 


{B'ij' + 2 Bi;") G 4 XX - ( -V'' - -g V'' - 2 ^^" ) G4X0 


\r B / ^ 

- 2BV'''G400 + ^ {5B'B^' + 6^2^-" - B'^|;' - 2B^") G 5 X + ) Gsx,^ 


2r2 

{B'll;' + 2 B^j") Gsxx + — + 45^^" + ^ + ) Gs,^ 

Zr^ r \ r r ' 


B 
2B2^'3 


G 


h(f)(f) 


(A3) 


c(*) _ 

^tt — 


G 2 x + G 30 + I (^4^ + + 2 ^;") Gax + ^ + S') G 4 X - + S^' + 25^") G 4 X 0 

(^SV' + 25^^" + ^V'') G4XX - ^ {SB'Bij' - BV - 2BV’" + 252,/^") Ggx - J 
^ + 2 SV 2 " + ^V'') G5X0 + (S>' + Gs XX , 


G 50 

(A4) 


£■1°) = 


2,/,/3 


G 2 + SV-'"G 2 x - B^'^Gs^ - ^ 


25^2 


BiIj^ 


/2 


4 2 / 

- + ^jG3x--(S- 


B-l 


O 4 - Bi/ 1 - + ^ ) G 44 , 


=>4 

r V A 

,/2 / Al 

3B^ 

r \ A 


2B- 1 

I- 

r 

3B- 1 


G 4 X 


BVq; + ^|G4X, 






A' , SV^A'^ 

— + - I G5X0--yGsxx , 


2r2 A 


(A5) 


£■(*) = 


— CT30 -I-;;-r^3X 


2 A 


BiP' A' 


2B A' 


-G. 


(3S - 1) Gsx + 


r A 


4X 


2BV^A' /4 A'\ 2BP;' 

Z<~r4(p(p - - -^CX4XX + ~ j '-^4X0- - — 


A' 1\ SVV ,1 

-Lr5X0 H— „ r, (j5xx + - 

A r ) 2 r^ r 


B-\ A' 

-B— 

r A 


G 54 , 


(A6) 
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= ujG 2 + [B'ij' + 2 B^") Ggx + f 


l'2B' A" B'A' y4'2_^2yl' 


A 


CJ 


B' 8 A' 

B r A 

B^i)' 


J + 2a;'' 


Ga- 


^Bib' + B'lh' + 2B^" + -Bib' ) a; - uj'Bib' 
A r ' 


G 


4(^ 


2A' , 4 5', A" , A'^ , B'A', A' „ 4 „ , 

7a^ + ;b^ ^ 


8ip' 25',, „ ,„ A',,. , 

—+ -g-^ +2i/i Iw 


G4x + B'^ip'^ ( ^ip' - ^ijj' + 2'0" - -if)' ) bj + ijj'w 

BA r 


— 2BiP'^u;Ga^^ + 


2„//3 


B^ljj 


-ijj'B' + i/j'B'^ + -BiP'iP" + 2 B^iIj" ) UJ - iijj'B' + 2Bil)") w 

r A r A ' 


G 4 X 4 . 

Lr4XX 


53 ^''^ 

4r 

B^iIj'^ 

4 

B^iP'^ 

4r 


5' 


A' 


6 


5 


,5'A 


A'- 


,A' 


A" 


- ( 5 -^ 1 /j' - —Ip' + 6ip" + -ip']uj'+{ 5 —— V;' - —pj' + &—tp” + 2—Ip' ) a; - 2iP'uj 


5 A 


A2 


A 


A 


G 5 X 


25 




,A' 2 , , 4 

/ I ./ / d' I 


i/;'5' + 2BiP" - iP' a;' - ip'B'— + 2B—iP" - -BiP'— + -iP'B' + -BiP" w 


A 


G. 


5X0 


{iP'B' + 2 BiP") oj' - — {ip'B' + 2BiP") UJ 


Gsxx + 


BiP' 


A' 8 \ 

— BiP' --BiP' - 4BiP" - 3B'iP' uj' 
A r 


A' BA''^ 9 A' A" fi S A' 

+ ( 3B'—iP' - —ip' + — -BiP' + 2B—iP' + -B'lP' + -BiP" + 4—BiP" ) uj - 2Bi\j' uj" 

AA^rA A r r A ' 


G 50 


2„/,/3 


B^iP 


2 A' 

-1- T 

r A 


UJ — UJ 


G 500 , 


(A7) 


= - ^^ujG3x - ujGs^ + 




2 A 


, ^^A" ^^A'^ ^ a; 


A' 25 

3—5V;' + — Ip' + B'lP' + 25V;" ) uj - uj'BiP' 
A r 


A +2S-r-25—+5' „ 
r A A A2 A 


Gax4> + 


- 5' 


^- 254 : 

r A 


- Buj" 


Gax 


BiP'{^iP' + 2iP"-iP'^']uj' 


A' 

--BiP' 


If I' 0 I" 2 / / \ 

~ +2ip - -Ip ) a; 


G 4 XX + 2ujGA(f><f> + 


B^PB'A'^, A' ^A" A'2 \ B^iP' „ 


iP'- 


AB' 

A^' 


+ 


4r 

B'^iP' {A'"^ , 5'A' , 

■—{-W* -Ta* 

B^iP'^ A' 

4r A 
5 

J 

5 

Y 

g^BA' /A' 


2r 


VB^ 

4r 

Gsx + 


6 V;' A' , 5' , , 

V-^A^ + ^ 5 ^ + 2 ^"j uj' 

_ ^^f _ 2 ^" + 2^' ) uj' 


A 


5 


A' ,,, 

6 A' 

2^r 

H-^ 

r A 

2A 



]a;'+( 


iP' 


-—iP' + -iP" uj 

r B r 


B'lP' A'lP' 


5 


2 A' A'2 A" 2 5' 5' A' 

r A +^A^"^T+r:B ■ 5 A + 


+ 2V’" j a; 

5' A' 

5 A 


G^xip 
G 5 XX 
a; 4- 2a; 


G 54 , 


A 4 A 


— o;'^ ^G5X0 — 2 G 4 XX H—^Gsxx^ , 


(AS) 
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c(0) _ 

— 


+ 


G 24 , + Btp' 

Bt/j' 

B 


B' 2 
2A^W r 
B' A' 4 2'iP' 

B A r ij;' 

tP" A' B' 2 


'll)'' 


+ H-h ^ 1 G 2 X + B-ip'^ G2x<j> — Bip 


Ip' 

G 30 - bV'" 


3A'B' 




/ ^/2 

~ ~A 
/ /1'2 

B^ip' 


A'B' 

2AB 


2v4' 

Ar 


2B' 

Br 


AAB 
B' 
Br 
2 ^ 


iB' 

Br 

A'B' 


B' 

T 

AV" 

Alp' 

2^)" 


AAB rip' 


BiP" 

Ip' 

AiP" 

'ip'r 

A'lP" 

2AiP' 


G 2 XX 
A" 


A'^ 

4A2 


2A' 

Ar 


Gsx 


Gsxx — Bip'^Gstjjtj, 


Br"^ J 


G, 


icf) 


A^r 


3 ,,3 /2A" 
(^47 




2A" 

5A'B' 

3A' 

3S' A' 

2A'0" 

20" , S' 

Ar 

ABr 

Aj-2 

Sr2 ' ABr'^ 

Aip'r 

ip'r"^ ' S 2 r 2 

A '2 

QA'B' 

3A' 

6 S' 8A'0" 

1 1 

80" 

S' 

20 " 

A 2 r 

ABr 

Ar 2 

Sr2 ' A0'r 

0 'r 2 

to 

to 

to 

Bip'r'- 


2 V’" 

Bip'r'^ 

Gixx 


Gax 


2„;y2 




A'^ 

742 ' 

2B' 


f ^- 

\ Ar Br 


2A'B' 

AB 

A'B' 


4A' 8 S' 

Ar Br 
A'tP" Af 


4 

■*" “2 
" 2 


SAV" 12V’' 


+ B^iP' 


2A" 

Ar 


/ ‘^ 4' 
S3V’/2 t 


B'^iP"^ 
B^'iP'^ 
B^'iP''^ 


V2Ar2 
A'2 


A'2 3A' 3A'B' 

A^r Ar^ ABr 

3A'2 ISA's' 


2AS A'V’' Ip'r 
3B' 

Br'^ 

A '2 

4ASr2 4A2Sr2 
SA'S' 7A'V’" 

L4V2 “ 2 Ar 2 “ 2 ASr 2 “ 77^ ^ 

A '2 ^// 7 ^/^/ 75 / 

~ 2ASr ~ 2Ar2 “ 2Br‘^ 

1 1V - 

Gs^^ S 0 - 


A' 


Sr2 A'V’' 

G4XX0 - SV'® 
2A'0" 20' 


Ip'r 

A'B' 


2„;y3 


G4X<^ + S 0' 


S' 2A'S0" 
^ + 


A' 4 

AT - 

A r 

2S0' 


G, 


Aip'r 0'r2 


4X00 


G 4 XXX 


ABr'^ 

A" 


S' 

A0'r ' ip'r'^ B'^r'^ 
3A'B' SA'iP" 


4A2r2 

A" 


2 ASr2 

A'B' 


AAB^r^ 

A'0" 


4AS2r2 2AB'ip'r' 


A' 


'^2 
5A'tP" 


Aip'r"^ 

Gsxx 

50" 


20 " 

tp'Br'^ 

A'lP" 

ABip'r"^ 


G 


50 


Gsx 


S' 


0 " 


2A2r 

A' 


2ASr2 
1 


Aip'r 0'r2 


G. 


5X00 


4.,y5 / 


A' 


S' 


+ 


A'0" 0" 


2ABr 2Ar2 2Sr2 A0'r ip'r"^ 


G5xx<p + S 0' 


A'BiP" 
\y4Ar2 ' 2Ar'^ip' 


A,„fA'B' 


ip' Br^ 


Gsxxx 


G^x4> 
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c(*) _ 

<-'A. - 
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A'BiP' 


G 


'34>4> 




2A 
B' 

Jb 

4y4'2 2A 


G 2 XX 


BA' 3BA'^ A'B' BA' 


-Bij}' 
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Ar 
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- H-h — 

r 2A ip' 

3A'B' A' 2A'^' 

Ar 


3x4, 
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Gsx 


A^r Ar ABr 
/2A'^B BA" A'B' 2 A'B 


2A ^ 

Ar 4A^ 
Aip'r ABr ^) 


A'^" 

2AiP' 


Gaxx 

4AB / 


I A^ 


A 

B' 

~B 


2A 

2^" 

Ip' 


Ar 


2B' 
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2(5-1) 


G. 


ix<t> 


^,,(2A B’ 2f" 4\ n2,/2/^6A' 2B' A'B' 
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; ,,3 fA'B' 2A'B'ijj" A'^B A'B\ 
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A'lP" 

Alp' 




Aip'r A^r Ar'^) 

2 ( 3A'2 A" 3A'B' 3A'2 A" A'B' \ 

V4A2r-2 “ 2Ar2 “ 4A^ ~ 4A^Br^ ^ 2ABr‘^ ^ 4AB^) 

/ A // -r^ A /O T~t n A I 7-»/ r» ^ / 7-» / // At 7-« / A ! ill 
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B 


\2Ar'^ 2A'^r‘^ 2Ar'^ 2A'ip'r'^ 

A" 5A'2 3A'B' A' 3B' 


4ABr^ 

3A' 


AiP'' 
2A'tp'r^ 
A'lP" 


A'^ 
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-B^iP 
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2jP 

~n—~2 
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^ \4ABr‘^ ^ 2AjP'r^ 4A^r^) ' 
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